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1. Introduction 

In the framework of non-abelian gauge theories at finite temperature, several effective descrip- 
tions have been pursued in order to overcome the infrared problems [[I]] connected with perturba- 
tive approaches to the fundamental theory. A rather successful technique is dimensional reduction 
[§, ||]. Thanks to the presence of different energy scales, induced by the finite temperature dy- 
namics of the original (3 + 1) -dimensional theory, an integration over the hard modes leads to a 
3d effective description which can then be solved in a non-perturbative way (e. g. by Monte Carlo 
integration). 

In the case of QCD, this technique loses its validity in the confined phase; however, one would 
want to devise effective methods to study the vicinity of the deconfinement transition: this is not 
a completely trivial task since the standard perturbative dimensional reduction does not retain the 
Z(N) symmetry of the original Yang-Mills theory [Q]. One can then follow a different strategy, 
namely writing down a general theory respecting the desired symmetry and then fixing the (many) 
couplings by matching with particular observables []|, g]. While for Si/ (2) the phase transition 
is captured correctly by such approaches [Q, 0], for the physically relevant SU (3) gauge theory a 
satisfactory fixing of all couplings is still an open issue. 

A different way to pin down a 3d effective theory is to employ lattice strong coupling expan- 
sions. This idea, first considered in [g], has been pursued by various authors @ [Tl], [jj, H HI 



and leads to theories with Polyakov loops as fundamental degrees of freedom. The contribution 
from spatial plaquettes was often neglected, a simplification which preserves the universal be- 



haviour of the theory; in Jig], instead, they were explicitly taken into account. Recent develop- 
ments including staggered fermions can be found in [|l7|]. 

The models proposed here systematically extend this approach by providing series for the 
effective couplings up to a certain order and are thus valid beyond the spatial strong coupling 
limit. As is to be expected from strong-coupling expansions, our results will have a finite radius 
of convergence, which is supposed to coincide with the deconfinement transition: in this sense, 
our effective formulation is complementary to weak coupling approaches. The effective actions we 
propose are subsequently studied by means of Monte Carlo integration, and the results are shown 
to lead to the correct order of the transition as well as to good estimates of the deconfinement point. 

2. Derivation of the effective theory 

2.1 General strategy and SU (2) 

Consider the partition function of a (3 + 1) -dimensional lattice gauge field theory at finite 
temperature ( T = J with gauge group SU (N) and Wilson's gauge action 



J [dU ] [dUi] 



exp 



4]>t/ p + trt/t) 



2N 



Finite temperature and the bosonic nature of the degrees of freedom imply the use of periodic 
boundary conditions in the time direction. 
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In order to arrive at an effective three-dimensional theory, we integrate out the spatial degrees 



of freedom and get schematically Q13Q 



[dU ] exp [-S eff ] ; 
—Ses = In j [dUi\exp 



XiSi +X2S2 + . . . . 



(2.2) 



We expand around j6 = and arrange the effective couplings X n = X n (f},N t ) in increasing order 
in j3 of their leading terms. Thus, the X n become less important the higher n. As we shall see, the 
interaction terms S n depend only on Polyakov loops 



Aft 



Lj = trWj = tr Y[U (xj,T) 



(2.3) 



With sufficiently accurate knowledge of the relations A n (j8,A^ T ), we are able to convert the cou- 
plings of the three-dimensional theory to those of the full theory. Determining the critical parame- 
ters X nx of the effective theory then gives a whole array of critical p c (N t ) for - in principle - all N x . 
In the following we calculate strong coupling, i.e. small j8, expansions of the leading X n . 

Since the calculations are quite similar for different numbers of colours, we now specialise 
our derivation to the simpler case of SU (2) and later provide the necessary changes for SU (3). For 
more details see [18]. Using the character expansion as described e.g. in [19, E9], the effective 



action according to Eq. (2.2) can be written as 



-S eff = ln f [dUi\H 

•J n 



r=iO 



(2.4) 



where the sum extends over all irreducible representations r with dimension d r and character % r . 



The expansion coefficients a r (j3) are accurately known [19] and in the following we use u = cif 
as expansion parameter instead of /3 for its better apparent convergence. The logarithm in this 



definition allows us to use the method of moments and cumulants Q21[], and we get the following 
cluster expansion 



-S eff = £ a(C)H^(xr,{Wj} 



III 



(2.5) 



*(xr,{Wj} 



[dUi]Y[dr p a rp Xr p {Up) , 

pex, 



where the combinatorial factor a(C) equals 1 for a single polymer X,- and —1 for two non-identical 
connected polymers. For clusters consisting of more than two polymers, a(C) depends on how 
these polymers are connected. Our task is then to group together all graphs yielding the same 
interaction terms up to some order in j8 , and this finally gives the strong coupling expansion of the 
corresponding effective coupling A„. 



3 



Effective Polyakov-loop theory 



J. Langelage, S. Lottini 




Figure 1: First graph with a nontrivial contribution after spatial integration for a lattice with temporal extent 
N z = 4. Four plaquettes in the fundamental representation lead to an interaction term involving two adjacent 
fundamental Polyakov loops L, and Lj, 



2.2 Leading order effective action 



The leading order result of the effective action has first been calculated in JT0| ] and corresponds 
to a sequence of N % plaquettes that wind around the lattice in temporal direction, cf. Fig. [jj. Its 
contribution is given by: 



<ij> 



(2.6) 



Hence, to leading order the first coupling of the effective theory is Xi(u,N T ) = u Nz . 

For additional terms of the series for X\ , we can use most of the graphs that also appear in the 
strong coupling expansion of the Polyakov loop susceptibility [^3|]. These corrections involve addi- 
tional plaquettes, are hence of higher order in u and we call these attached plaquettes decorations. 
Carrying out the calculations, we get the following results through order u 10 in the corrections 
relative to the leading order graph: 



Ai(m,2) = w 2 exp 
Ai(m,3) = w 3 exp 
Ai(m,4) = w 4 exp 
X]_(u,N T > 5) = M*exp 



, , 4 „ fi 134 o 49044 l0 

2 4u 4 - Su 6 + u u 

3 405 



3 l 4 u 4 -4u 6 + —u s 



X\Au A -Au 6 + —u 



36044 
405 ' 

37664 
405 



.10 



10 



. , 4 , fi 140 o 36044 
A^ T ( 4u 4 - 4u 6 + — — — u - -jjft-u 



10 



(2.7) 



For smaller A^ T some graphs do not contribute since the temporal extent of their decoration is > N t 
so that they do not fit into the lattice. 



2.3 Higher order terms 

There occur several types of higher order graphs: larger numbers of loops involved, Polyakov 
loops at distances larger than one and Polyakov loops in higher dimensional representations. We 
begin by considering powers of the leading order term. Inspection of higher order terms shows that 
one can arrange a subclass of these terms in the following manner 

£ fw r ^ + & L 3_ Jj = £ m( i +Al L iLy ) . (2.8) 

<ij> V z / <ij> 



4 



Effective Polyakov-loop theory 



J. Langelage, S. Lottini 



To see this, one calculates the corresponding graphs with LfLj or Z|Ly, and the combinatorial 



factor a{C) of Eq. ( |2.6| ) gives the correct prefactors for the series to represent a logarithm. 

Next, let us consider couplings pertaining to next-to-nearest neighbour interactions. These 
appear once additional plaquettes are taken into account. Naively, the leading contribution should 
correspond to a planar graph with Polyakov loops at distance two. However, this graph is precisely 
cancelled by the contribution of the nearest-neighbour graph squared and its associated combina- 
torial factor —1. The leading non-zero contribution therefore comes from L-shaped graphs and is 
given by 

h(u,N t )S 2 =N t (N t - 1) M 2 ^ +2 £L*L, , (2.9) 

M 

where we have two additional spatial plaquettes and we sum over all pairs of loops with a diagonal 
distance of \ 2a, abbreviated by [Jd]. With the same steps leading to Eq. (2.8), we finally arrive at 
the SU (2) partition function 

Z= ( [dW] J] [1 + MULj] n [1 + hL k Li] . (2. 10) 

J <'j> [kl] 

Finally, we include some remarks about the Polyakov loops in higher dimensional representations. 
Consider, e.g., the adjoint Polyakov loop: the leading order term emerging from a strong coupling 
expansion is 

KS a = V N ^ £ Xa(Wi)Xa(Wj) , V = \? + \vt + ^~U 6 + . . . 

<ij> " 

and hence X a ~ u 1Nt , which is formally of lower order than the coupling %z. To next-to-leading 
order (valid for all A^ T > 2) we have 

X a = ^(l+N^ + ..y (2.11) 



Effects of higher representations have also been investigated in the literature [ 14, O, 25]. 



2.4 The effective action for St/ (3) 

In the case of SU (3) the same steps as for SU (2) apply. The only difference we have to keep 
in mind is that SU (3) also has an anti-fundamental representation and consequently there is also 
a complex conjugate Polyakov loop variable L*. Thus we get the one-coupling and two-coupling 
partition functions 

Zi = / [dW] n [y + h{UL)+L*Lj)] , (2.12) 

J <ij> 

Z 2 = [ [dW] n [l+Ai {L i L*+L*L j )]Y\[\+X 2 (L k L* l +L* k Li)] . (2.13) 

J <U> [kl] 
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The effective coupling X\{u,N x ) is obtained as (for this gauge group we consider only even values 
ofiV T ): 

Ai (2, u) = u 2 exp [2 (4m 4 + 12m 5 - 18m 6 - 36m 7 



219 o 1791 n 830517 I0 



Ai (4, u) = u 4 exp [4 (4m 4 + 12m 5 - Uu 6 - 36m 7 



295 o 1851 n 1035317 , n 

H m s H u H u 

2 10 5120 

Ai (# T > 6, b) = m Wi exp [# T (4m 4 + 12m 5 - 14m 6 - 36m 7 

295 o 1851 n 1055797 m 

H m H m H M 

2 10 5120 

For the first terms of the next-to-nearest neighbour coupling X2(N t ,u) we find 

A2(2,m) = m 4 2m 2 + 6m 4 + 31m 6 , 
A 2 (4,m) = m 8 [i2m 2 + 26m 4 + 364m 6 ] , 

A2(6,m) = M' 

A 2 (/V T > 8,m) = u 2N * [N Z (N Z - 1)m 2 ] 
while the leading coupling of adjoint loops is (valid for Af T > 2) 



(2.14) 



.12 



30m 2 + 66m 4 



A fl = v^(l+A^ + ...y 



9 2 9 3 81 4 

V = — M M H M + . 

8 8 32 



(2.15) 



(2.16) 



3. Numerical simulation of the effective theories 

3.1 The one coupling model 

For the purpose of numerical simulations, a further simplification is achieved by using the 
trace of the Polyakov loops for the path integral measure as degrees of freedom (complex numbers, 
\L X \ < 3, instead of matrices), and rewrite the one-coupling partition function for SU (3), Eq. (2. 12), 

Z = (J] j dL,V 5rff ; Seff = - £ log(l + 2A 1 ReL i L}) -£V V . (3.1) 

x J <ij> x 

The potential term V x is the Jacobian induced by the Haar measure of the original group in- 
tegration; rotating the matrices to the diagonal form diag(e' e ,e~'( 9+ ^), with |0|, |0| < k, we 
have p6|]: 



V x = ^log(27 - 18IL,! 2 + 8ReZ^ - |L X | 4 ) . 
The integration measure actually used in our simulation then takes the form 



dL x e 



+ 71 



dd x e v < . 



When working on the SU (2) theory, — 2 < L x < +2 is a real number and we simply have 

"+2 



£ *L x e v *, V x = l -\og{A-L 2 x ). 



(3.2) 



(3.3) 



(3.4) 
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N 3 = 06 
N, = 08 
N s = 10 



0.03 
0.025 

0.02 
0.015 

0.01 
0.005 




N. = 06 ' — »— , ' 
N s = 10 J 
- N s = 12 d * T 




/ ,4 

»* r/ 1 

... ' 


\ K 







0.0925 0.093 0.0935 0.094 0.0945 0.095 0.0955 



Figure 2: Left: Distribution of L for small and large X\ on a lattice with N s = 6 and M = 1. Middle, Right: 
Expectation value of |L| and its susceptibility. The vertical line marks the infinite-volume transition. 



3.2 A "sign problem" and its solution 

Our numerical approach will be a straightforward Metropolis local update algorithm; however, 
the Boltzmann weights to consider are in the form exp(log(l + 2AiRe(L,L*))): for high enough 
couplings, they can be also negative, thus spoiling the update technique (the partition function 
being, overall, still positive). In the SU (2) case, the threshold coupling X T = 1/4 is well beyond the 
phase transition, so that in practice there is no problem around criticality, but the SU (3) threshold 



of 1/9 is very close to the transition and a direct numerical investigation of the model as in Eq. ( p7T[ ) 
is impossible. 

Our approach to overcome this problem is the following: we Taylor-expand the logarithm in 
the effective action to some order M in powers of q = AiReL,//- (undoing the resummation as in 
Eq. (fT8|)), obtaining models free of the problem: 

32 



AM) 
5 eff 



I 

<ij> 



2q — 2q + —i 



(3.5) 



In this way we can identify a critical point for each M and look for their convergence as M — > oo. 
Also, we can compare to the SU (2) case where the M = oo value is directly calculable. 

3.3 Phase structure, critical coupling and finite size analysis 

Our first task is to establish the phase structure of the effective theory, where we focus on the 
physically interesting case of SU (3). Based on the global Z(3) symmetry of the model, one expects 
spontaneous breaking of that symmetry for some critical value of the coupling X\, c . Fig. |2] shows 
the behaviour of the field variable L as a function of X\ . As expected from the 4d parent theory, 
there is indeed a transition from a disordered or mixed phase, with values of L scattering about zero, 
to an ordered phase at large coupling where the three Z(3)-phases are populated separately. In the 
thermodynamic limit, one of these vacua will be chosen and the symmetry is broken spontaneously, 
(L) = for X\ < Ai f and (L) ^ for X\ > X\ jC . Correspondingly, the expectation value of \L\ rises 
abruptly at some critical coupling Xi jC , as shown in Fig. || (middle). On a finite size lattice, the 
phase transition is smoothed out, non-analyticities are approached gradually with growing volume, 
as the figure illustrates. 

The critical coupling, X\ x , is located via finite-size scaling. After identifying a pseudo-critical 
&i,c(Ns) for a number of finite systems, the relation 



-l/v 



(3.6) 
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data i 
first-order fit 
TD limit 0.094238(10) 
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Figure 3: Left: Position of the minimum of the Binder cumulantB(£') for SU(3),M = 1, for different lattice 



sizes. The horizontal line is the thermodynamic limit resulting from the fit to Eq. (3.6). Right: Behaviour 
of B m i n (N s ), along with its thermodynamic limit obtained with the ffiN^) scaling law and the independent 
estimate Boo from the \L\ histogram. Also the second-order limit value 2/3 is shown. 



is used, with v = 1/3 for the SU (3) first-order transitions and, in the SU (2) case, the 3d Ising value 



V = 0.63002 [27]. Numerically, we found satisfactory results with data produced in just a few days 
on a desktop PC. 

For the definition of the pseudo-critical coupling, one can look at the energy E = — S e ff/Ai 
(neglecting the potential term) or derived quantities, but in general, due to the nonlinearity of 
5 e ff in the coupling, we preferred to look at the average modulus |L|; one can then define Ai c as 
the minimum/maximum of the associated Binder cumulant/susceptibility, which indeed featured a 
more robust scaling: 



S(|L|) = 1 ~3lt# ; X(\L\) = {{\L\-(\L\)) )■ (3-7) 

3.4 Critical coupling and order of the transition for SU (3) 

The truncated theories with M = 1,3,5 were simulated on lattices with spatial sizes N s = 
6,8, 10 (plus N s = 12 for the M = 1 theory). For each volume, ~ 30 values of the couplings are 
sampled by ~ 10 6 update sweeps each. Measurements were taken every ~ 30 updates. 

Regardless of the truncation order M, the SU (3) theories display a first-order transition; among 
the associated features, we found very long thermalisation times <=c expfcA^ 3 ) as is expected for 
tunnelling phenomena (Fig. |]): for instance, a system with size N s = 16 would require, around 
criticality, ~ 10 6 update sweeps to thermalise. 

First we consider the model with M = 1 . The first-order nature of the transition is established 
by fitting the pseudo-critical couplings to the scaling law, Eq. (3i)), with v = 1/3, see Fig. || 



(left). The behaviour of the minimum B m [ n of B(|L|) is a further confirmation; this quantity, as 
demonstrated in [|^, ^9h, scales as B min (N s ) = + B^N~ 3 + ff(N~ 6 ), with a thermodynamic 



limit which is smaller than the second-order value 2/3, 

!-—( J£Jj_ _ 

3 12V|L| 2 |Lh 



2 1 /|L|i |L| 2 \ 2 
Boo = -- — (^-TT-r) , (3-8) 
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2 p 
1.5 - 




Figure 4: Left: Behaviour of L with Monte Carlo time for two N s = 6 trajectories in the SU(3) M = 1 theory 
with X\ = 0.0935. Right: Histogram for L, obtained from 60 such trajectories. The tunnelling between the 
central and the three broken-symmetry vacua is apparent. 



data from B( L ) 
data from v{ L ) 
M=infinity, B( L} 
M=infinity, x( L ) 



1 2 3 4 5 12345678 

M (truncation order) M (truncation order) 

Figure 5: Truncation-dependence of the critical points. Left: in the 51/(3) case, the data points re- 
fer to three system sizes, while the lines mark the extrapolated critical points; the latter are found at 
A Lc = 0.094238(10), 0.10635(11), 0.10403(28) for M = 1,3,5 respectively. Right: for SU(2), pseudo- 
critical points for N s = 8 obtained from B(|L|) and #(|L|) for a variety of truncations (data points) and 
compared with the untruncated M = °° values (lines). 



N s = 6 data i — i — i 
N s = 8 data > « 
N s =10data . 
Critical. M=1 

Critical, M=3 

Critical, M=5 


0.24 
0.23 
















£ 0.22 








0.21 
















0.2 



with |L|i and \L\i the two local maxima of the \L\ double-peaked histogram. A direct comparison 
between the results for Boo from scaling analysis and from the location of \L\i shows an agreement 
within two standard deviations, the residual discrepancy being probably due to neglecting higher- 
order N~ 6 corrections. 

In the next step we need to investigate the behaviour of the models with higher M. Again 
we observe first order transitions, which become sharper with increasing M. Moreover, finite- 
size effects are stronger for higher M, Fig. |5| (left). The critical couplings identified for the M = 
1,3,5 effective theories in the thermodynamic limit are also quoted there. Judging from these 
three values, the series seems to be rapidly converging, with only ~ 3% difference between M = 
3,5. The residual difference between this estimate and the M = oo critical coupling is completely 
subdominant compared to the other systematic errors contributing to the final results. Also, the 
direct comparison with the SU (2) case below, where the M = oo data are directly available, supports 
a rapid convergence, Fig. || (right). 
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3.5 Critical coupling and order of the transition for SU (2) 

In this family of theories the transition is second-order; with much less relaxation problems 
(e.g. 4000 steps for N s = 16), larger lattices (up to N s = 28) were available. With the same ap- 
proach as for SU(3), the nature of the transition was confirmed by: (a) X\ tC (N s ) scaling with the 
3d Ising critical index, (b) Binder cumulant analysis approaching 2/3 for large systems, and (c) 
\L\ histogram inspection, where a single peak continuously moves to the right as the coupling is 
increased. All inspected values of M yielded the same features. Moreover, here a direct compari- 
son with the M = oo untruncated model is possible, and shows that a rapid convergence is indeed 
realised (Fig. ||, right); in particular, we found 

X hc (M = 1) = 0.195374(42) ; X hc (M = oo) = 0.21423(70) , (3.9) 

which indicates quite small systematic deviations due to choosing one particular truncation. 

3.6 Two-coupling models for SU (3) 

In this section we study the influence of including a second coupling. We consider two possi- 
bilities: the first one is switching on the interaction between next-to-nearest neighbours. The SU (3) 



version of Eq. (2.10) reads: 



z = ( n / dL n c 1 + 2 ^i Rez ^p nc 1 + 2^^;)^ . (3.10) 

* J <ij> [kl\ 

We remark that now there are two terms suffering from the above-mentioned sign problem: a 
truncated expansion is then needed in both, and the two truncation parameters {M\,M2) should be 
chosen in a consistent way, for all N t , with respect to the power in u we want to keep. We adopted 
the choice (3, 1) after checking numerically that higher values of M2 give negligible differences in 
the results. 

In the other model, we allow the nearest neighbours to interact also in the adjoint representation 
as described before. The partition function in this case (with the adjoint part already truncated at 
M2 = 1) is given by 

II / dL *) EI (1 +2AiReLiL}) [I e^'^CaW^eExVi , (3.11) 

X <ij> <ij> 

with the adjoint trace Tr^W = | Tr W | 2 — 1. Also in this case, the truncation (3,1) was employed. 

In these two-dimensional parameter spaces, there is a critical line separating the symmetric 
and the broken phases. However, for a given N T , only a one-dimensional manifold in this space 
represents the image of the original gauge theory, since both couplings are functions of the sole u. 
The strategy was then to identify the shape of the critical line and find, for each temporal lattice 
extent, the intersection with the curve enforcing that particular value of N T . 

In both models, the critical lines were found by interpolation after locating 1 1 critical points 
at as many fixed values of the second coupling; it turned out that a linear parametrisation was good 
enough in describing them (within our precision, finite-size effects were practically invisible): 

X hc = a + bX 2 with a= 0.10628( 8), b = -1.891( 4) . (3.12) 
X lc = a + bX a with a= 0.10637(15), b = -1.422(22) . (3.13) 



10 



Effective Polyakov-loop theory 



J. Langelage, S. Lottini 




Figure 6: Critical line in the two-coupling space, determined from Dashed lines give the parameter 

space representing a4d theory with fixed A^ T . Left: (hi, fa). Right: (fa, fa). 





M= 1 


M = oo 


4d YM 


3 


2.15537(89) 


2.1929(13) 


2.1768(30) 


4 


2.28700(55) 


2.3102(08) 


2.2991(02) 


5 


2.36758(40) 


2.3847(06) 


2.3726(45) 


6 


2.41629(32) 


2.4297(05) 


2.4265(30) 


8 


2.47419(22) 


2.4836(03) 


2.5104(02) 


12 


2.52821(14) 


2.5341(02) 


2.6355(10) 


16 


2.55390(10) 


2.5582(02) 


2.7310(20) 



Table 1: Critical couplings j3 c for 51/(2) from two effective theories compared to simulations of the 4d 



The value of a was always, as expected, compatible with the estimate for the critical point of the 
M = 3 one-coupling theory. 

By plotting these critical lines and the family of curves coming from requiring a given N T , 
one sees that the latter accumulate towards vanishing second-coupling as N t increases (Fig. ^): 
this implies that the effect of including those interactions is less and less important at finer lattice 
spacings: only at very low values of temporal extent does the inclusion of a second coupling make 
any visible difference. 



4. Mapping back to 4d Yang-Mills 

Having established the critical couplings for our effective theories and tested their reliability, 



we are now ready to map them back to the original thermal Yang-Mills theories by using Eqs. ( |2.7| , 
114]). In Tables |, | we collect the values for the critical gauge couplings, j8 c , obtained in this way 
from the effective theories and compare them to the values obtained from simulations of the full 4d 
theories for SU (2),SU (3), respectively. 

The agreement is remarkable in all cases, with the relative error of the effective theory re- 
sults compared to the full ones shown in Fig. [7| The comparison of alternative truncations of the 
logarithm shows once more that it has almost no influence on the accuracy of the final result, as 
described earlier. Interestingly, there appears to be a 'region of best agreement' , with the deviation 
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N t 


M= 1 


M = 3 


M 1 ,M 2 (A 2 ) = 3,1 


M u M 2 {l a ) =3,1 


4d YM 


4 


5.768 


5.830 


5.813 


5.773 


5.6925(002) 


6 


6.139 


6.173 


6.172 


6.164 


5.8941(005) 


8 


6.300 


6.324 


6.324 


6.322 


6.0010(250) 


10 


6.390 


6.408 


6.408 


6.408 


6.1600(070) 


12 


6.448 


6.462 


6.462 


6.462 


6.2680(120) 


14 


6.488 


6.500 


6.500 


6.500 


6.3830(100) 


16 


6.517 


6.528 


6.528 


6.528 


6.4500(500) 



Table 2: Critical couplings j3 c for St/ (3) from different effective theories compared to simulations of the 4d 



M-1 results > 
M.infinity results * 



M=3 ■ 
(M,,M z (y . 3,1 



Figure 7: Relative error of j3 c predicted by the effective theories when compared to simulations of the 4d 
theories, for 51/(2) (left) and 51/(3) (right). 



growing both for small and large N r . We ascribe this to the fact that there are two competing sys- 
tematic errors, as discussed earlier: the validity of the strong coupling series for a given coupling 
A; is better the smaller j8 and hence N T , whereas the truncation of the next-to-nearest neighbour 
interactions gains validity with growing N t . In particular in the case of SU (3), there appears to 
be a cancellation of the two kinds of systematics, rendering the effective description better for the 
original theory on finer lattices. 

The strong-coupling series was inspected both by comparing the resulting j8 c from series of 
different depth and by Pade analysis, and we observe a satisfactory convergence. It was also found 
that the error due to the truncation of the strong-coupling series is much larger than that from 
neglecting higher couplings. 

One can also compare the results presented here with those from the inverse Monte Carlo 
approach, where the effective theory is found in a completely non-perturbative way; inspection of 
the SU (2) case |]T4[], in particular, shows that the abrupt change of curvature in the inverse Monte 
Carlo function X\ (j8 ) at the critical point is not captured by our strong-coupling approach, which is 
consistent with j8 c marking the radius of convergence also for the series expansion of the effective 
coupling X\ . Thus, the inverse Monte Carlo approach has a wider range of validity whereas the 
series approach furnishes analytically known mappings between the full and effective theories. 
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5. Conclusions 

We have derived, by means of strong coupling expansions, an effective description for lattice 
pure gauge theories at finite temperature which respects explicitly the requirement of centre sym- 
metry and has only scalar Polyakov loop variables as degrees of freedom. Moreover, due to the 
dimensional reduction involved, the ^-dependence is encoded completely in the maps from the ef- 
fective to the original coupling j3 , whose expansion can be extended, in principle, to higher orders. 
We have also considered interaction terms other than the leading one, namely a next-to-nearest 
neighbour interaction and an adjoint-representation coupling term. 

Our Monte Carlo approach to the models, while requiring modest computational resources, 
confirms the expected nature of the symmetry-breaking transition for both SU (2)- and SU (3)-based 
effective formulations, and allowed us to predict the critical point j3 c of the original 4d thermal 
gauge theories with an accuracy within a few percent for a variety of values of N r . Particular 
attention was devoted to estimating the effect of employing different approximations, with quite 
stable answers in support of the good convergence of the series, and of neglecting higher-order 
interaction terms, which again does not have a strong effect on the final answers especially at finer 
lattices. 

An extension of the present work could be the study of SU (N) gauge theories with N > 3 



(cf. Q34| ] and references therein), which can be performed much in the same way as the cases 
examined here; even more intriguing is the possibility to keep the theory simple while getting a 
step closer to physical QCD, i. e. by introducing fermions and finite baryon density, for instance by 



employing a hopping parameter expansion [23, 35] 
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